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COMBINATION  EQUATIONS  USED  TO  CALCULATE  EVAPORATION 
AND  POTENTIAL  EVAPORATION^ 

J.  F.  Bartholic,  L.  N.  Namken,  and  C.  L.  Wiegand^ 
INTRODUCTION 


Some  of  the  best  known  and  most  recent  evaporation 
equations  were  derived  from  basic  energy  balance  and 
turbulent  transfer  principles  using  consistent  nomen- 
clature. Since  few  of  the  equations  differentiate  between 
soil  evaporation  and  plant  evaporation,  the  flux  of  water 
vapor  away  from  any  surface  will  simply  be  referred  to 
as  evaporation. 

In  figure  1  all  arrows  (fluxes  of  energy)  toward  the 
canopy  are  positive,  while  energy  fluxes  away  from  the 
canopy  are  negative;  that  is,  /?„  is  positive  while  E,  H, 
and  G  are  negative  values  in  this  figure  (all  symbols  are 
defined  in  the  appendix). 


Figure  1.- Schematic  diagram  of  fluxes  of  energy  to  and  from 
the  canopy. ^^^_^_^__^_^B^. 
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DEVELOPMENT  OF  EQUATIONS 


ENERGY  BUDGET  EQUATIONS 

The  vertical  energy  balance  for  a  cropped  surface  can 
be  given  by 


Rn+E  +  G  +  H  =  Q. 


(1 


This  equation  is  depicted  in  figure  1  where  net 
radiation  (Rn)  is  adding  energy  to  the  canopy  and  it  is 
being  repartitioned  into  soil  heat  flux  (G),  evaporation 
(£),  and  sensible  heat  (//).  With  well-watered  fields  and 
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moderate  climates,  the  ratio  of  evaporation  to  net 
radiation  (£"//?„)  is  nearly  equal  to  one  while  with  a  dry 
soil  surface  or  stressed  plants  EIR„  <1  is  expected  with 
the  ratio  frequently  as  low  as  0.5.  If  sensible  heat  is 
advected  into  a  plot,  the  ratio  can  be  greater  than  one 
and  reported  values  as  high  as  1 .5  are  not  uncommon. 
The  following  assumptions  are  made  in  equation  1 : 

(1)  The  surface  is  large  and  uniform,  that  is,  there  is 
no  horizontal  divergence. 

(2)  Energy  of  photosynthesis  is  negligible.  It  is 
usually  only  1  to  2  percent  of  net  radiation,  although 
Lemon  has  stated  it  can  go  as  high  as  5  percent  during 
the  most  active  growth  stage  of  corn  (7).^ 

(3)  Energy  stored  in  the  crop  volume  can  be 
neglected.  This  is  generally  valid  (2^),  except  possibly 
during  a  short  period  near  sunrise  and  sunset. 


Underscored  numbers  in  parentheses  refer  to  Literature 


Cited,  p.  11. 


In  assumption  number  1 ,  it  is  difficult  to  put 
numbers  on  what  is  a  sufficiently  large  uniform  surface 
to  provide  negligible  horizontal  divergences.  In  1958, 
Elliott  (4)  analyzed  height  and  fetch  requirements,  and 
reported  that  the  maximum  height  (z,)  for  suitable 
measurement  with  a  95  percent  profile  adjustment  could 
be  found  approximately  by  the  equation  z,-  =  0.75A'°'' 
for  A'  as  a  fetch  distance  in  centimeters.  This  equation 
worked  if  the  fetch  was  greater  than  10  meters  (m.). 
Rider  (14)  evaluated  temperature,  moisture,  and  profile 
variations  downwind  from  a  tarmac  surface  for  the  first 
16  m.  over  a  grass  surface.  The  measurements  showed 
that  the  greatest  variation  occurs  in  the  first  10  to  1  5  m. 
after  a  discontinuity  is  encountered. 

Dyer  and  Crawford,  however,  reported  data  for  a 
well-watered  site  surrounded  by  dry  fields  in  which  they 
were  able  to  measure  horizontal  variations  as  far  as  230 
m.  downwind  over  the  well-watered  field  (3).  They  also 
showed  how  evaporation  and  sensible  heat  transfer 
would  vary  from  the  leading  edge  inwards  distances  of 
150  to  200  m.  Generally,  most  of  the  measurable 
horizontal  variations  near  the  surface  have  dissipated  in 
the  first  100  m.  Instruments  above  the  crop  need 
roughly  200  m.  of  fetch  for  each  meter  above  the  crop 
surface  for  a  90  percent  profile  adjustment. 

R„  is  generally  measured  directly,  but  it  also  can  be 
obtained  from 

Rn    =    Qs    -    Qsr  ^  Qa-     Qar  "     Qbs  I^] 

if  each  of  the  components  of  /?„  are  known. 

TRANSFER  EQUATIONS 

The  flux  of  heat  through  a  solid  (//j)  can  be 
determined  by  the  equation 

where  K  is  the  coefficient  of  heat  conduction  (thermal 
conductivity)  for  the  solid.  This  coefficient  remains 
relatively  constant  over  a  wide  temperature  range  for 
solids. 

The  turbulent  vertical  transfer  of  heat  in  air  is  given 
by 


H  =  pc,K,  [If  +  r] 


[3] 


where  T  is  the  dry  adiabatic  lapse  rate  dT/dz  =—0.01 
C/m.  and  can,  therefore,  be  neglected  in  many  applica- 
tions of  [3],  since  dT/dz  over  the  earth's  surface  is 
generally  10  to  100  times  greater  than  P,  except  at 
sunrise  and  sunset.  The  temperature  in  equation  3, 
strictly  speaking,  should  have  a  bar  over  it,  since  it  is  a 


time  averaged  value.  However,  since  the  temperatures 
and  water  vapor  measurements  referred  to  in  this  paper 
are  all  time  averaged,  no  bars  are  used. 

The  equation  for  the  turbulent  flux  density  of  latent 
heat  is  commonly  seen  in  the  following  two  forms  for 
gradients  of  specific  humidity  and  vapor  pressure, 
respectively. 


E  =  LypK, 


dz 


^      ^     0.62    de  ,41 


where  0.62  is  the  ratio  of  the  molecular  weight  of  water 
to  air.  The  flux  of  momentum  or  Reynolds'  stress  can  be 
given  as 


T  =  pK, 


dz 


pu^' 


[S] 


where  the  stress  is  assumed  constant  with  height. 

If  the  coefficients  for  turbulent  heat,  vapor,  or 
momentum  transfer  in  the  air  {K^^  K^,  and  K„^, 
respectively)  were  as  constant  as  that  for  heat  conduc- 
tivity in  a  solid,  the  calculation  of  evaporation  or 
sensible  heat  flux  or  both  would  be  simple.  However,  the 
transfer  coefficient  above  a  surface  depends  on  wind- 
speed,  surface  roughness,  and  air  stability. 

The  roughness  affects  the  size  of  the  "turbulent 
parcels"  and  a  hot  surface  affects  the  buoyancy  and 
shape  of  the  parcels.  Because  the  transfer  mechanism 
above  the  soil-plant  interface  is  such  a  dynamic  and 
complex  process,  the  main  problem  concerning  equa- 
tions 3,  4,  and  5  has  involved  the  manner  in  which  the 
transfer  coefficients  vary.  In  developing  many  of  the 
following  evaporation  equations,  the  manner  in  which 
the  transfer  coefficients  are  handled  is  of  key  impor- 
tance. 


EQUATIONS  FOR  EVAPORATION 

Bo  wen  Ratio  Method 

Bowen  (J_)  used  a  combination  of  the  energy  budget 
[  1  ]  and  the  transfer  equations  3  and  4.  In  this  method 
the  ratio  of  H/E  is  used  where 


B  -- 


H       CpP    Kh  dTldz 


0.62  L^K^  de/dz 


Assuming  Kfj  =  Ky^  and  that  the  gradients  are  measured 
over  the  same  interval,  then  upon  integration 


B  =  H/E  =  7  (Ar/Ae) 


[6] 


where  7  =  Cp  Pj Q.62 
rearranging  gives 


Ly.  Substituting  [6]   in  [1]   and 


-  JRn  +  g) 

1   +  T    (AT/Ae)' 


[7] 


or 


" = m 


ph   iea  -  e^) 


Asterisks  are  used  beside  the  equation  number  of 
commonly  used  evaporation  equations. 

Equation  7  generally  works  well  when  the  evaporative 
flux  is  large  relative  to  sensible  heat  flux.  Under  these 
conditions  AT  is  small  and  Ae  is  large,  so  errors  in  the 
second  term  in  the  denominator  have  little  effect  on  the 
calculated  value  for  E.  This  evaporation  equation  has 
been  used  extensively  by  Frischen  and  van  Bavel  (5). 
Pruitt  (^),  Tanner  (22),  and  others  and  is  generally 
adequate  for  determining  evaporation. 


+   /  P^    {(e'a  -  e^)  -  (^   -  e^)}       [12] 


substituting  (e^  -  e^,)  from  [9]  in  the  first  term  on  the 
right  side  of  [  1 2]  gives 


H=     r|j    E+-^    ph    \(e',   -e^)-(e'a-eo)\        [13] 


Vapor  Pressure  Deficit  Combination  Formula 

Slatyer  and  Mcllroy  (18^)  and  Penman  (1^)  derived  an 
equation  for  evaporation  that  used  [1],  [3],  and  [4]. 
Equation  3  for  this  derivation  is  written  in  the  integral 
form 

H  =  pCphiTa  -   Ta)  [8] 


where  h  =  Ki,/z  and  [4]  is  written  as 


Z,vP  0.62 
E  =- n h  (ea  -  eg). 


[9] 


Symbols  with  a  subscript  "o"  refer  to  surface  values. 

The  slope  of  the  saturation  vapor  pressure  curve  for 
water  (A)  (19)  can  be  written  as  A  =  9e  /bT  or  relative  to 
the  terminology  being  used,  the  slope  is  closely  approxi- 
mated by 


A  = 


ie'a  -  e'o) 


[10] 


Substituting  [13]  into  [  1  ]  yields 


A 


A  +  7 


{- 


(R    +    G) 


p-^h     \i'4  -ea)-{e'o  -  ^o)\\ 


141 


This  is  a  general  form  for  the  combination  method 
involving  vapor  pressure  deficits.  In  [14]  (e^  —  e^)  is  the 
water  vapor  saturation  deficit  at  some  height  in  the  air 
while  {c'q  —  Cq)  is  the  saturation  deficit  at  the 
evaporating  surface. 

With  different  surface  moisture  conditions,  [14]  can 
be  modified.  If  the  surface  is  wet,  then  no  saturation 
deficit  exists  and  e'^  =  e^,.  Using  this  restraint  in  [14] 
yields  the  following  equation  for  potential  evaporation 


£-p  = 


A  +  7 


\-  (R„  +  G)-p^h{e'a-e^)^   _    [15]* 


for  small  differences. 

Rearranging  [10]  and  substituting  it  in  [8]  gives 


^  =  P^    M<  -e'o) 


A 


11 


where    the    primed    values    refer    to    saturated    vapor 
pressure  at  the  corresponding  subscripted  temperature. 

Upon  substituting  the  mathematical  equalities  e'^  = 
^a  +  (^'a  ~  ^a)  and  e'^  =  e^,  +  (ej,  -  e^),  we  can  write 
[11]  as 


-H 


ph  (ea   +  e'a   -  ea  -  eg   -  Co     +    ^o). 


In  practice  T^  is  seldom  known,  so  in  using  equations 
14  and  15  an  additional  assumption  is  frequently  made. 
It  is  assumed  that  T^  can  be  used  for  Tq  and  A  then  is 
given  by  the  slope  of  the  saturated  vapor  pressure  curve 
for  water  at  T^.  Ta  is  also  used  in  the  calculation  of  7. 
This  assumption  is  beUeved  to  cause  little  error  in  the 
calculation  of  evaporation  for  well-watered  surfaces 
(24). 

A  further  modification  of  this  equation  is  possible  if 
evaporation  is  occurring  from  a  wet  surface  with  infinite 
fetch,  since  in  this  case  the  air  near  the  surface  is 
saturated  and  (e^  -  e^)  =  0  and  [15]  reduces  to 


A 
^P  A  +  7 


(Rn  +  G). 


[16] 


Equation  16  can  be  put  into  the  form  of  the  Bowen 
ratio  equation  7,  where  E  becomes  Ep  from  a  wet 
surface  with  infinite  fetch: 


Ep 


+  G) 


-(Rn     +     G) 


1 

A 


1    +  7 


(Tg     -     TqY 

(e'a    -    e'o) 


171 


Evaluation  of  the  Transfer  Coefficient  h 

In  evaluating  the  transfer  coefficient  h  the  relation- 
ship is  really  being  determined  between  h  and  wind- 
speed,  surface  characteristics,  and  stability  of  the  air. 
The  development  of  this  relationship  generally  starts  by 
assuming  an  equality  between  the  transfer  of  momentum 
and  that  of  water  vapor. 

With  this  basic  assumption  in  mind,  a  study  of 
information  available  relating  to  the  wind  profile  is 
helpful.  The  slope  of  the  wind  profile  varies  inversely 
with  height  above  the  ground.  The  relation  is  expressed 
by 


du 
dz 


[181 


where  the  constant  of  proportionaUty  is  the  frictional 
velocity  (u^)  divided  by  von  Karman's  constant  (k^ 
0.4).  Integration  of  [  1  8]  from  u  =  o  at  z^  to  u^  at  z-d 


unstable 


neutral 


Figure  2. -Schematic  representation  of  windspeed,  distribution 
in  and  above  a  vegetative  system  with  different  stability  regimes. 
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where  z^  is  the  roughness  length  and  d  is  the  zero  plane 
displacement  as  shown  in  figure  2.  Values  for  Zq  and  d 
can  be  obtained  from  wind  profile  data  using  a  least 
squares  fit  of  the  data  to  a  wind  profile  equation  such  as 
[19]  or  by  graphically  fitting  the  best  curve  to  the  data. 
The  graphical  fitting  of  the  data  is  shown  in  figure  3. 
Windspeeds  above  a  140  centimeter  (cm.)  tall  cotton 
canopy  with  neutral  temperature  conditions  are  shown 
in  part  A  and  the  same  data  is  shown  in  part  B  with  a  d 
value  of  60  cm.  The  height  is  plotted  on  a  logarithmic 
scale  for  both  A  and  B.  Extrapolation  of  the  profile  to  u 
=  o  gives  a  z^  value  of  1 9  cm. 
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Figure  3. -Wind  profile  data  above  a  140-cm.high  cotton  canopy: 
A ,  under  neutral  conditions;  and  B,  with  a  d  value  of  60  cm. 


Unfortunately  [19]  holds  only  for  neutral 
(adiabatic)  conditions  as  shown  in  figure  2  and  not  for 
unstable  (diabatic,  lapse)  conditions  or  during  stable 
(subadiabatic.  inversion)  conditions.  Thus,  the  log  law 
holds  only  for  neutral  conditions  and  does  not  ade- 
quately incorporate  other  stabiUties  which  it  must 
include  to  be  of  general  applicability  (25). 

Equation  18  can  be  generalized  by  including  a 
stability  function  as  follows: 


dz 


[201 


Upon  integration 


0/z  3z 


21 


One  form  given  for  0  that  is  frequently  referenced 
was  suggested  by  Monin  and  Obukhov  (9)  where 


1  +  a  -  . 


[22] 


and/.  =-TpCpU^  KkgH) 


in  which  T  is  an  absolute  temperature  and  a  is  an 
empirical  constant.  Another  form  for  0  involving 
Richardson's  number  (Ri)  for  modification  of  the 
logarithmic  wind  equation  with  diabatic  conditions  is 


=  (1   +  bRi)", 


and 


n-    -    S      ,  iz    -    d)  -2 

Ri   -  —   \n     (Ta-To)  Ua        z. 


[23] 


[24] 


Lettau  (8^)  has  tabulated  the  different  values  for  a  and  b 
obtained  in  six  research  studies.  Other  forms  of  0  are 
numerous  and  Sellers  (j_6,  p.  153)  states,  "At  least  eight 
different  'diabatic  profiles'  have  appeared  in  the  litera- 
ture, each  one  supposedly  fitting  the  observations." 

The  values  for  a  and  b  that  are  most  likely  to  be 
recommended  are  a  =  -0.25  and  b  =-18.  Thus  [23]  is 


=  (1   -   ISRO' 


[25] 


Equation   25    has  been  derived  independently  by  five 
authors  and  has  been  termed  the  "KEYPS"  profile. 

Since  the  transfer  coefficient  /j  is  in  [9] ,  which  is  [4] 
integrated,  we  will  continue  the  evaluation  of  h  by 
reintegrating  [4] 


'/; 


z-d 

E  I       i//Cw  az  =  / 


J   eo 


be 


which    upon    integration    with    /    =    .62    L^pIP,    and 
rearranging 


o  ~  ^o)\J  '^       UKy,   dz\ 


E  =  i(e 


but  K^  as  a  function  of  z  is  not  known. 


However,  replacing  r—  in  [5]  with  [20]  gives 
oz 


K,y,   =  ku. 


[26] 


[27] 


Using  the  assumption  stated  above  that  A^^  -  K^^  allows 
the  substitution  of  [27]  in  [26] 


E  =  jku*  (e^  -  e 


-[/::■•■]■■ 


[28] 


Since   u,    may    be    difficult    to  evaluate,  it   can  be 
substituted  from  [21  ]  which  gives 


E  =  k'^  Uajiea  -  Cq) 


[/ 


z-d 


-,  -2 


bz 


[29] 


Setting  [9]  equal  to  [29]  gives 


/ 


z-d  -"2 

-     OZ 


[30] 


A    form   of  h  found  in  the  literature  (J_6)  can   be 
obtained  by  simply  letting 


0  ^  0  -   1    ^  I 
z  z  z 


in  [30]  which  upon  integration  gives 


k'    u„ 


-d  ^2 


h'"^] 


where 


rz-d 

J     Zn 


-     1 


3z, 


[31] 


[32] 


The  expression  for  h  in  [31]  can  be  simplified  for 
neutral  stability.  With  neutral  stability,  0  =  1  and  [31] 
red'ices  to  the  logarithmic  law  for  the  wind  profile. 
Likewise,  [20]  reduces  to  [  18] ,  and  [25]  reduces  to  (p  = 
I  since  Ri  reduces  to  zero  as  shown  by  [24]  because  in 


£•=£■„    + 


pCr 


neutral  conditions  T^, 
rewritten  as 


Tg.  With  *  =  0,  [31]  can  be 


k'    u„ 


[331 


where  Ca  is  the  drag  coefficient. 

Over  water,  well-watered  vegetation,  or  wet  soil,  the 
temperature  profile  is  frequently  near  adiabatic.  Under 
these  conditions  and  when  working  near  the  crop 
surface,  h  as  evaluated  from  [33]  is  generally  believed  to 
work  satisfactorily. 

Under  stable  conditions,  [22]  has  frequently  been 
found  to  best  represent  the  atmospheric  stability  cor- 
rection factor;  however,  since  stable  conditions  exist 
most  frequently  at  night  when  evaporation  is  small  or 
negligible,  stable  conditions  can  be  considered  minor. 

An  exception  to  this  general  statement  may  exist 
when  oasis  (advective)  conditions  exist  and  hot  air  blows 
over  a  cool,  well-watered,  rapidly  evaporating  crop. 

For  unstable  conditions,  /;  is  evaluated  from  [31] 
which  requires  the  integration  of  [32].  Using  numerical 
integration,  Lettau  (S)  developed  the  table  of  h  values 
for  different  stability  conditions,  and  Tanner  (22)  has 
also  published  the  same  table.  Fuchs  and  Tanner  (6) 
reported  that  /;  values  calculated  with  [31]  under 
unstable  conditions  above  a  drying  soil  averaged  25 
percent  larger  than  /;  values  obtained  from  [33] . 


A  Generalized  Combination  Method  by  Fuchs  and 
Tanner 

Fuchs  and  Tanner  have  proposed  a  combination 
formula  relating  E  to  Ep  in  which  Ep  and  then  E  are 
evaluated  (6).  Their  equation  can  be  obtained  by  using 
[15]  in  [14]  and  using  the  identity  (e'^  -  Cg)  =  {e'o  ~ 
ej)~(^o  -«'a)  to  yield 


E     = 


E     +      -^ 

"  (A +  7) 


[    (^o    -   ^a)     ~     (^o    ~  ^a)  ] 


Then  substituting  from  [9]  for  —h  (e^  ~ e^)  gives 


[34] 


A  +  7 


[m^;,  -e. 


)    +    El 


L^p  .62 


Rearranging  gives 


[A  +  7]  pCp 


[35] 


36] 


The  authors  used  [31  ]  to  obtain  /;  in  [36] .  G  and  Tq 
are  the  only  surface  measurements  required  in  this 
equation.  The  larger  Tg  is,  the  larger  the  second  term  on 
the  right  becomes  and  the  smaller  E  becomes,  remember- 
ing that  evaporation  always  has  a  negative  value  when 
the  vapor  flux  is  away  from  the  surface.  Physically  this  is 
to  be  expected  since  as  E  decreases,  with  the  same  /?„, 
the  energy  must  be  partitioned  differently.  An  increasing 
Tq  allows  this  to  occur  since  H  and  G  away  from  the 
surface  can  increase  as  well  as  Qiy^.  Fuchs  and  Tanner 
found  a  good  comparison  between  the  E  values  from 
[36]  and  those  from  the  Bowen  ratio  method  for  a  soil 
with  a  nearly  dry  surface. 

For  smooth  soil  or  water  surfaces,  the  surface 
temperature  can  generally  be  used  as  Tq.  With  crops  or 
rough  surfaces,  however,  it  should  be  kept  in  mind  that 
To  is  the  air  temperature  at  c?  +  z^,  since  this  is  the  lower 
limit  used  in  the  integration  of  the  equations  used  in 
deriving  most  of  the  evaporation  equations.  Thus,  a  true 
To  or  eg  can  be  obtained  only  by  an  extrapolation  of  air 
temperature  or  water  vapor  pressure  profiles  to  the 
height  {d  +  Zo).  Thus,  equation  36  works  best  for 
nonvegetative  surfaces. 


Evaporation  Equations  Using  Transport  Resistance 
Concepts 

Several  investigators  (1_0,  20)  consider  the  resistance 
concept  as  shown  in  figure  4  of  value. 

Evaporation  equations  for  the  fluxes  involved  can  be 
written  as 


E  =  / 


E  =  j 


E  =  / 


To height  a 


ei T± vapor  source 


Figure  4. -Resistive  network  model  of  the  transport  of  heat  and  vapor  in  the  atmospheric  surface  layer  and  through  the  substrate 

surface  layer. 


Here  r^  =  l/h  for  the  air  resistance.  However,  if  a 
canopy  is  involved,  an  additional  resistance  for  the 
canopy  should  be  included.  Also,  a  resistance  value  for  a 
laminar  flow  surface  layer  and  leaf  stomata  for  cropped 
surfaces  may  need  to  be  included  in  some  cases. 
Although  this  approach  may  be  appealing,  there  are 
possible  inaccuracies  in  the  model  since  evaporation 
from  the  soil  as  well  as  the  canopy  takes  place  from  a 
vertically  distributed  array  of  sources  as  opposed  to  the 
one  horizontal  vapor  source  postulated  in  resistance 
models.  Tanner  and  Fuchs  (23)  have  discussed  some  of 
the  problems  with  resistance  models. 

EQUATIONS  FOR  POTENTIAL  EVAPORATION 

Evaporation  equations  provide  a  measure  of  water 
being  lost  with  a  given  atmospheric  demand.  Potential 
evaporation  equations,  on  the  other  hand,  yield  evapora- 
tion that  would  occur  with  a  given  atmospheric  demand 
from  a  "wet  surface";  that  is  a  water  surface  or  a  field 
recently  watered  where  there  is  no  restriction  on  the 
availability  of  water  for  evaporation.  In  the  equations 
discussed  in  this  paper,  if  the  development  of  the 
equation  involves  having  e^  be  the  saturated  vapor 
pressure  at  T^,  the  equation  assumes  a  wet  surface  and 
the  equation  predicts  potential  evaporation. 

The  Combined  Concept 

One  potential  evaporation  equation,  [15] ,  was  devel- 
oped in  the  discussion  of  vapor  pressure  deficit  formula. 
Van  Bavel  (24)  used  this  equation  and  evaluated  h  with 
equation  33,  which  gives 


A 

^P   '    A  +  7 


\     -iRn     +     G) 


With  this  concept  no  measurements  of  surface  tempera- 
ture are  required;  however,  this  may  cause  some  error  in 
A.  Possible  effects  of  this  error  in  A  were  evaluated  and 
found  to  cause  only  a  very  small  error  in  Ep. 

The  use  of  [37]  by  other  authors  has  cast  some 
doubt  on  the  vahdity  of  the  form  of  /;  used  in  this 
equation  during  windy  conditions.  In  particular  with  a 
value  of  Zq  of  1  cm.  or  greater,  the  second  term  becomes 
very  important.  Skidmore  (17)  reported  that  in  Kansas 
for  average  daily  windspeeds  of  0.88  and  2.26  m/sec, 
measured  at  45-cm.  height,  the  wind-dominant  term 
(second  term)  contributed  33  and  113%,  respectively,  as 
much  as  the  radiation-dominant  term  to  the  total 
calculated  potential  evaporation.  For  these  two  consecu- 
tive days,  the  ratio  of  Ep  from  [37]  to  net  radiation  was 
0.98  and  1.60,  respectively,  while  the  ratio  of  Ep  over  E 
from  the  Bowen  ratio  were  1.16  and  1.62. 

Parmele  (U)  found  that  during  midday  equation  37 
tended  to  overestimate  evaporation  from  bare  wet  soil 
by  10  to  25  percent  and  evaporation  was  also  generally 
overestimated  for  perennial  ryegrass.  However,  these 
were  truly  overestimates  only  if  the  measured  evapora- 
tion was  occurring  at  a  potential  rate.  Parmele  found 
when  using  24-hour  average  values  in  [37]  that  calcu- 
lated and  measured  daily  evaporation  were  very  close. 
Using  Zq  values  suggested  by  van  Bavel  on  crops  like 
well-watered  alfalfa  and  daily  values  of  the  required 
parameters,  Jensen"*  found  total  calculated  evaporation 
for  the  season  to  be  about  two  times  net  radiation, 
which  was  not  realistic.  Using  a  value  of  0.5  to  1.0  cm. 
for  Zq,  however,  yielded  values  that  compared  favorably 
with  measured  values  with  alfalfa.  Good  results  from  the 
equation  have  been  reported  for  both  hourly  and  daily 
values  over  open  water,  wet  soil,  and  alfalfa  when 
measurements  taken  at  2  m.  were  used  (24). 


[in    ^1'    (4  -  ea)\.      [37] 


Personal  communication  from  Marvin  E.  Jensen. 


Potential  Evaporation  Using  Surface  Temperature 

Conaway   and   van    Bavel   {2)    proposed   the  use  of 
equation  9  with  h  determined  from   [33],  which  gives 


E^n   =  /^^    ^a 


[■"^] 


{ea-eo).       [381* 


This  method  is  based  entirely  on  mass  or  vapor  transfer 
equations    and    does    not    involve    the    energy    balance 


equation.  In  this  approach  surface  temperature  is  re- 
quired, and  Conaway  and  van  Bavel  suggest  using  infrared 
thermometer  techniques  for  determining  its  value.  They 
define  e^  as  surface  vapor  pressure;  however,  it  probably 
would  be  more  proper  to  define  e^  in  this  potential 
evaporation  equation  as  the  saturated  vapor  pressure  at 

To- 

This  equation  gave  accurate  values  of  £"_  from  a  wet 
soil  surface  both  on  an  hourly  and  daily  basis. 


SUMMARY 


In  deciding  which  evaporation  equation  will  work 
best  in  a  particular  project,  several  criteria  should  be 
evaluated,  including  the  following: 

(1)  Is  evaporation  or  potential  evaporation  to  be 
evaluated? 

(2)  Are  instruments  required  to  make  the  measure- 
ments available,  and  can  all  the  variables  required  in  the 
model  be  measured  with  sufficient  accuracy''  Can  data 
reduction  and  calculations  be  handled  without  delay  and 
considerable  man-hours? 

(3)  Is  the  evaporation  equation  based  on  a  physical 
or  empirical  model  and  has  the  equation  been  suffi- 
ciently evaluated?  If  not,  are  lysimeters  or  other  means 
available  to  adequately  check  the  equation?  Only  equa- 
tions derived  basically  from  physical  models  have  been 
discussed  in  this  publication;  however,  Rosenberg  has 
recently  reviewed  the  broad  scope  of  evaporation  equ- 
ations (IJ). 

(4)  What  assumptions  and  boundary  conditions  have 
been  used  in  developing  the  equation,  and  can  these 
conditions  be  met  in  this  particular  research  project? 

Certainly  other  criteria  could  be  suggested,  but  these 
four  should  be  sufficient  for  choosing  the  best  evapora- 
tion equation  to  be  used  in  a  particular  project. 

Whether  an  evaporation  or  potential  evaporation 
equation  is  needed  will  depend  entirely  upon  the  type  of 
research  project.  Certainly  a  water  surface  as  well  as 
many  well-watered  crops  and  soils  will  evaporate  at  or 
near  their  potential  rates.  Evaporation  will  be  less  than 
potential  when  any  restriction  on  the  availabiUty  of 
water  for  evaporation  exists.  Frequently  the  Ep  equa- 
tions are  less  complex  than  actual  evaporation  equations. 

All  of  the  required  measurements  for  the  evaporation 
equations  derived  can  be  made;  however,  the  ease  with 
which  the  different  measurements  can  be  made  varies 
considerably.  Of  the  energy  flux  values  required,  net 
radiation  (/?„)  is  probably  the  easiest  to  measure,  but 
good  instrument  calibration,  maintenance,  and  periodic 
recahbration  are  required.  Since  jR„  is  large,  small 
percentage  errors  in  it  greatly  affect  E.  Soil  heat  flux  G 

n 

has  a  much  smaller  value  than  /?„  (G  <T;r.  generally)  so 


larger  percentage  errors  in  G  can  be  tolerated.  G  can  be 
obtained  from  soil  heat  flux  plates,  which  have  been 
calibrated  in  the  soil  with  a  moisture  content  near  that 
expected  during  the  experiment,  or  by  calculation  from 
soil  temperature  and  moisture  measurements  and  a 
knowledge  of  some  physical  properties  of  the  soil. 

Vapor  pressure  of  the  air  (e^)  is  frequently  obtained 
from  wet-  and  dry-bulb  measurements,  LiCl  probes, 
infrared  adsorption,  or  dew  cells,  and  any  of  these 
methods  can  be  sufficiently  accurate.  All  methods  for 
determining  water  vapor,  however,  require  extreme  care 
and  considerable  maintenance.  Saturated  vapor  pressure 
for  a  smooth  wet  surface  e'^  is  obtained  with  ease  if 
surface  temperature  is  known,  since  e'g  =  f(Tg).  If  £■  < 
Ep  then  e'^  "^^  e^,  and  e^,  is  very  difficult  to  determine 
accurately. 

Air  temperature  (Tg)  is  relatively  easy  to  measure  but 
To  may  be  more  difficult  because  direct  contact 
measurements  of  surface  temperatures  are  seldom 
acceptable  under  natural  conditions.  Remote  deter- 
mination via  infrared  thermometry,  however,  can  be 
used  for  smooth  surfaces  to  determine  T^  if  emissivity 
of  the  surface  and  background  radiation  in  the  appro- 
priate wavelength  interval  are  known.  Infrared  deter- 
mination of  T(,  makes  possible  the  use  of  equation  36  to 
determine  E  and  equation  38  to  evaluate  £"„. 

In  the  Bowen  ratio  method,  gradients  in  the  air  of  Ae 
and  Ar  may  be  required.  Over  the  intervals  for  which 
the  gradients  are  determined-usually  an  interval  of  only 
50  to  possibly  200  cm. -the  differences  in  temperature 
and  vapor  pressure  may  be  quite  small;  thus,  con- 
siderable accuracy  is  required.  With  a  well-watered 
surface  Ae  may  be  large,  while  AT  will  usually  be  small 
and  should  be  known  to  at  least  ±0.02  C.  For  a  surface 
or  crop  that  is  quite  dry,  AT"  will  generally  be  large, 
while  Ae  will  be  small.  Since  Ae  appears  in  the 
denominator  of  equation  7,  when  it  is  a  small  fraction  it 
greatly  affects  the  calculated  value  for  E.  Because  of  the 
accuracy  required  in  obtaining  these  gradients,  thermo- 
piles or  other  very  sensitive  instruments  for  obtaining 
temperature  differences  are  used,  and  the  sensors  at  the 
two  levels  are  periodically  interchanged  to  eliminate 
instrument  bias. 
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Evaporation  equations  that  require  an  evaluation  of 
the  turbulent  transfer  coefficient  h  require  a  continuous 
measure  of  u^  and  periodic  wind  profile  measurements 
with  at  least  six  sensitive  cup  anemometers.  Wind  profile 
measurements  under  neutral  conditions  are  required  for 
the  evaluation  of  d  and  z^.  Also,  for  some  crops  Zg 
appears  to  be  a  function  of  windspeed,  so  profiles  with  a 
variety  of  wind  conditions  may  be  required  as  well  as 
periodic  measurements  as  the  surface  changes  in  rough- 
ness or  as  the  crop  grows.  Since  a  number  of  values  have 
been  reported  for  z^  and  d  in  the  literature  for  various 
crops,  use  of  "book"  values  for  Zq  and  d  may  be 
possible. 

For  rapid  data  reduction  and  analysis,  a  complete 
system   analysis  incorporating  the  entire  system  from 


initial  fabrication  or  purchase  of  sensors  through  experi- 
mental setup,  data  recording  systems,  equipment  main- 
tenance, computer  computations,  printout,  and  plotting 
is  a  must. 

Various  empirical  evaporation  equations  have  been 
used  extensively  and  some  have  been  very  valuable  in 
providing  estimates  of  evaporation.  Nearly  without 
exception  the  empirical  equations  do  best  in  locations 
similar  in  chmate  to  where  the  equations  were 
developed;  in  other  cUmates,  these  equations  may  yield 
unrealistic  estimates  of  evaporation.  Evaporation  equa- 
tions based  on  physical  models  should  be  generally 
applicable  to  all  areas,  if  the  models  are  sound  and  the 
conditions  used  in  developing  the  equation  are  met. 
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TABLE  OF  TERMS  AND  UNITS 


APPENDIX 

Symbol 


Symbol 

Cn 


Description 


d 
E 


ei 


G 

8 
H 


["■^^] 


drag  coefficient 


(dimensionless) 


specific   heat  of  air  at  constant  pressure  (cal 

g.-'  °c.-') 


height  displacement  parameter  (cm.) 

latent    heat    flux    density    (evaporation)    (cal. 
cm."^  min."' ) 

potential   evaporation   flux   density  (cal.  cm."^ 
min."' ) 

water  vapor  pressure  (mb) 

vapor  pressure  in  the  air  at  height  a  (mb) 

saturated  water  vapor  pressure  at  Ta  (mb) 
water  vapor  pressure  at  T,-  (mb) 
water  vapor  pressure  at  surface  (mb) 
saturated  water  vapor  pressure  at  Tg  (mb) 
soil  heat  flux  density  (cal.  cm.'^  min."* ) 
acceleration  of  gravity  (980  cm.  sec."^) 
sensible  heat  flux  density  (cal.  cm."^  min."' ) 


K 

Kh 

k 
P 

Qa 

Qar 

QbS 

Qs 

Qsr 


Description 

transport  coefficient  for  heat  and  water  vapor 
(cm.  min."' ) 

=  0.62  L^pip  (cal.mb"'    cm."^) 

thermal  conductivity  (cal.  min.''  cm."'  °C."' ) 

transfer  coefficient  for  momentum  (cm) 

transfer  coefficient  for  heat  (cm.^  min."' ) 

transfer    coefficient    for    water    vapor    (cm.^ 
min."' ) 

von  Karman  constant  (^0.4)  (dimensionless) 

latent  heat  of  vaporization  (cal./g.) 

atmospheric  pressure  (mb) 

incoming  thermal  radiation  flux  density  from 
the  atmosphere  (cal.  cm."^  min."' ) 

reflected    thermal    radiation   flux   density  (cal. 
cm."^  min."' ) 

emitted    thermal    radiation    flux    density    (cal. 
cm."^  min."' ) 

incoming  solar  radiation  flux  density  (cal.  cm."^ 
min."' ) 


reflected  solar  radiation  flux  density  (cal.  cm 
min.~' ) 


-2 
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Symbol  Description 

q  specific  humidity  (dimensionless) 

R,j         net  radiation  flux  density  (cal.  cm."'^  min."') 

r^  transport    resistance    between  the  surface  and 

height  a  {r^  =  Ijh)  (sec.  cm."' ) 

r;,  internal  resistance  of  the  surface  layer  to  the 

transport  of  heat  (sec.  cm."' ) 

r,-  internal  resistance  of  the  surface  layer  to  the 

transport  of  water  vapor  (sec.  cm."' ) 

T  temperature  (  C) 

Ta  air  temperature  at  height  a  (  C) 

Tq  surface  temperature  (°C) 

u  windspeed  (cm./min.) 

u^  windspeed  at  height  a  (cm./min.) 

u  frictional  velocity  (cm./min.) 

2  height  in  the  air  from  the  soil  surface  (cm.) 

2/  height  above  the  crop  to  which  the  profile  has 

been   modified  to  95   percent   of  its  final 
value  (cm.) 

Zq  roughness  height  (cm.) 

a  constant  in  Monin  and  Obukhov  wind  profile 

equation  (dimensionless) 


Symbol 


Description 
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4> 


CnP 

(mb  °C."' )   psychometric  constant 


Z,v   .62 
r  dry  adiabatic  lapse  rate  (  C.  cm."' ) 

'z-d 


J""  111 

J       Zr^ 


dz  (dimensionless) 


0  stability  profile  influence  function  (dimension- 

less) 

p  density  of  air  (g.  cm.    ) 

A  slope  of  the  saturation  water  vapor  curve  versus 

temperature  (mb  °C.'' ) 

A^         vapor  pressure  difference  between  two  levels  in 
the  air  (mb) 


AT        temperature  difference  between  two  levels  in 
the  air  (°C.) 

T  flux  density  of  momentum  (g.  cm."^  min."' ) 


LIST  OF  MAJOR  EQUATIONS 

Energy  balance  equation 


Rn  +  E  +  G  +  H    =    0 


Turbulent  heat  transfer  equation 
H  =  pCpKh 


■■  [^ '  '] 


or 


H  =  pCph   (Ta   -  To) 


Latent  heat  transfer  equation 

dq 
E  -  Ly^pK^   -       -  L^pK^  ■ 

.62    de 
P     9z 

or 


Ly,p  .62     ,     ,  . 


Bowen  ratio  method 
-  {Rn   +  G) 


E  = 


1    +  7 


Ar 

Ae 


Vapor  pressure  deficit  combination  formula 

A 


A 


+  7       I 


{Rn    +   G) 


[1] 


[3] 


[8] 


[4] 


[9] 


[7] 


C 


[14] 


P-^    >^     \(<   -ea)-  ie'o   "  ^o)  jl 


A  generahzed  combination  method  by  Fuchs  and  Tanner 


E  = 


^-■'  -  ^p  , ,,.  -  ,^, 


A 


'^'      A 


:36i 


where 


h   =  k^  Ua        *    +     \n 


[31 
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and  Ep  is  calculated  t  om  the  vapor  pressure  deficit  where 

equation  15. 


2        r      (z-^)   T     _ 
The  combined  concept  (for   -p)  ^  ~  ^    ""     \^"       7^        J        "  <^a  "a  [33] 

A      r 
Ep    =— —{R„  +  G)  Potential  evaporation  using  surface  temperature 


pCp 


2'a  -e,)J  [151*  ^P   =  i'^'^^a    [l"^   1 


-2 


h     (4  -  ej  [151*  Ep   =  jk    Ua       1"——  iea   -  e^)  [3%] 
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